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Abstract 

We present exact results, at next-to-leading order in renormalization- group im- 
proved perturbation theory, for the Wilson coefficients appearing at order 1/uiq 
in the heavy-quark expansion of heavy-light current operators. To this end, we 
complete the calculation of the corresponding two-loop anomalous dimension ma- 
trix. Our results are important for determinations of \V u b\ using exclusive and 
inclusive semileptonic B decays. They are also relevant to computations of the 
decay constant /b based on a heavy-quark expansion. 
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1 Introduction 



Hadronic matrix elements of flavor- changing currents are of paramount importance to 
the phenomenology of weak decays. They appear, e.g., in the theoretical description 
of semileptonic B decays. These processes give direct access to the elements \V c b\ and 
\V u b\ of the quark-mixing matrix, provided the hadronic matrix elements are known 
from some nonperturbative computation. For inclusive semileptonic decays, the matrix 
elements can be calculated using a short- distance expansion in the large 6-quark mass. 
Improving the theoretical accuracy in the determinations of \V u b\ is currently one of the 
most important challenges in heavy-flavor theory. 

Relevant to the determination of \V u b\ are hadronic matrix elements of the weak 
current W7 a (l — 75)6. In general, local current operators composed of a heavy quark Q 
and a light quark q exhibit an interesting behavior under renormalization at scales below 
the heavy-quark mass. Large logarithms of the type a s ln(mg//x) arise from the exchange 
of gluons that are hard with respect to the light quark but soft with respect to the heavy 
quark. Since such gluons see the heavy quark as a static color source, their effects 
can be investigated systematically in the framework of the heavy-quark effective theory 
(HQET), which provides a systematic expansion of matrix elements in inverse powers and 
logarithms of mq In the HQET, the four-component heavy-quark field is replaced 
by a velocity-dependent two-component field h v satisfying j)h v = h v , where v is the four- 
velocity of the hadron containing the heavy quark. Operators in the effective theory 
have a different evolution than in usual QCD. For instance, whereas the vector current 
q~f a Q is conserved in QCD (i.e., its anomalous dimension vanishes), the corresponding 
current q^ a h v in the HQET has a nontrivial anomalous dimension, which governs its 
renormalization-group (RG) evolution for scales below tuq. 

As an example, consider the heavy-quark expansion of the vector current (neglecting 
the light-quark mass m q ) [Q, [J 

1 10 

- qi a Q - £ CiQi) J. + — ^Bjifj) Qj + 0(l/m 2 Q ) , (1) 

•=1,2 Zm Q j=l 

where the symbol = is used to indicate that this is a relation that holds after taking 
matrix elements on both sides. (In a regularization scheme with anticommuting 75, the 
expansion of the axial-vector current is obtained by replacing q — > —(775 in the HQET 
operators.) J x = qj a h v and J 2 = qv a h v are the effective operators entering at leading 
order in the expansion. The operators Qj appearing at next-to-leading order (NLO) 
form a basis of dimension-4 operators in the HQET that closes under renormalization. 
The standard choice of these operators is (using the notation iD = id + g s A and (iD)i = 
-id + g s A) 

Qi = q^iiph 
Q 2 = qv a HPK 
Q 3 = qiD a h v 



Q A = q(ivDyj a h v , 
Q 5 = q(ivDy V a h v , 

Q 6 = q(iD a yh v , (2) 
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and 



Qi 
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i J d 4 xT{J li2 (0), h v {iDfh v (x)} 



d 4 xT{J li2 (0),|^a^G^^(x)} 



(3) 



In addition to the local operators Qi,.s, there appear four bilocal operators Q7...10 built 
from time-ordered products of the leading-order currents with an insertion of the "kinetic 
operator" h v (iD) 2 h v or the "chromo-magnetic operator" ^h v a tl yG liy h v . These operators 
appear at order I/tuq in the effective Lagrangian of the HQET. Part of their effect is to 
take into account the expansion of the external hadron states in terms of states of the 
effective theory jl]. The Wilson coefficients of these nonlocal operators are given by the 
products of the coefficients of their component operators, i.e., 



£7,8 (aO = Ci.aO) , B 9t io(fi) = Ci, 2 (a*) CmagO) 



(4) 



where C mag (/^) is the coefficient of the chromo-magnetic operator in the HQET La- 
grangian, which has been calculated at NLO in 0]. The kinetic operator is not renor- 
malized ||. 

The calculation of the Wilson coefficients Cj(/i) and Bj(fi) at NLO in RG- improved 
perturbation theory requires one-loop matching at the scale \x = mq and two-loop 
anomalous dimensions in the effective theory. The leading-order coefficients Ci(fi) have 
been computed long ago || [7]. The calculation at order 1/toq is, however, much more 
complicated. The leading logarithmic corrections were obtained in f2|. First steps to- 
wards a NLO analysis were made in ||, where the one-loop matching calculations were 
computed, and general arguments based on the reparameterization invariance |5j of the 
HQET were used to derive the exact relations (valid to all orders in perturbation theory) 



bm = cm , bm = I bm = cm 



(5) 



Moreover, it was shown that the 10 x 10 anomalous dimension matrix governing the 
mixing of the operators Qj under RG transformations has the texture 



where 



7 



^7hi 7a \ 

7b 
V 7c 7£, / 



(6) 



7hi = diag(7i,7i,7i) , 

7o = diag(7i, 71, 7i + 7ma g , 7i + 7 

{lB)ij = 71 $ij + $i3 (lA)3j ■ 



(7) 
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Here 71 and 7mag are the anomalous dimensions of the leading-order heavy-light currents 
J12 and of the chromo-magnetic operator, respectively. The matrix 



1a = 



/ -2( 72 + 73) 2 73 2 72 \ 



V -(72 + 73) 73 72 / 



(8) 



describing the mixing of the local operators Q1...3 i n t° Q4...6 has been calculated at two- 
loop order in ||. Explicit expression for the two- loop anomalous dimensions will be 
presented in Section |5|. 

It was observed in || that there are numerically large discrepancies between the one- 
loop and leading-logarithmic approximations to the Wilson coefficients Bj(fi), which 
could only be resolved at NLO in RG-improved perturbation theory. The missing piece 
in the NLO analysis is the calculation of the anomalous dimension matrix -yc describing 
the mixing of the bilocal operators Q7...10 into the local operators Q4...6. In this paper, 
we report the calculation of this matrix at two-loop order. It requires the evaluation 
of two-loop tensor integrals in the HQET, which are infrared (IR) singular when one of 
the external lines is taken on-shell. A general algorithm for computing such integrals 
has been developed in For our analysis we will generalize the above discussion and 
consider an arbitrary Dirac structure V of the current qTh v in the time-ordered products. 
At the end we will study the particular case of vector currents in detail. 

The reader interested only in the final results of our calculation but not its technical 
details can continue with Section ||, where we present the final expressions for the Wilson 
coefficients Bj(n) at NLO and briefly discuss some phenomenological applications. 

2 Operator mixing 

Our task is to calculate, at two-loop order, the mixing of bilocal operators, consisting 
of time-ordered products of a local current q Th v with either the kinetic operator or the 
chromo-magnetic operator, into local dimension-4 current operators with a derivative 
acting on the light-quark field. Here V can be an arbitrary Dirac matrix, depending 
on the Lorentz structure of the original current that is expanded in terms of HQET 
operators. 

We start by constructing, in the HQET, a basis of the relevant operators that mix 
under renormalization. Since ultimately our interest is in the matrix elements of these 
operators between physical hadron states, it is sufficient to consider gauge-invariant 
operators that do not vanish by the equations of motions. For the case where the bilocal 
operator contains the kinetic operator, such a basis consists of the two operators 



O 



kin 



q(ivDyTh. 




(9) 
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If the bilocal operator contains the magnetic operator, the basis consists of the three 
operators 

2 mas = -\ q (iDjvyJ r (1 + i) a^K , 

0™ ag = ijd 4 x T{gIX(0), | ^(^(a;)} . (10) 

For each of the two cases, we define a matrix Z of renormalization constants, which 
absorb the ultraviolet (UV) divergences in the matrix elements of the bare operators, 
by the relation Oi = J2j %ij Oj.barc- The matrix 7 of the anomalous dimensions, which 
govern the scale dependence of the renormalized operators, is given by 

in) 

dm /1 

In a minimal subtraction scheme, the renormalization constants are defined to remove 
the 1/e poles arising in the calculation of the bare Green functions with insertions of 
the operators Oi in dimensional regularization, i.e. in d = 4 — 2e space-time dimensions. 
Hence, 

z = i + Y,-z {k) . (12) 

k=Q 6 

The requirement that the anomalous dimension matrix be finite in the limit e — > 
implies the relations |§ 

n dZ^ <9Z (2) dZ^ / m /3(a s )\ 

7 = 2u s —- , a s — = a s - Z« + , (13) 

oa s oa s oa s \ a s J 

where P(a s ) = da s /dln/x 2 is the (3 function. The first equation shows that the anomalous 
dimension matrix can be obtained from the coefficient of the 1/e pole in Z, whereas the 
second one implies a nontrivial constraint on the coefficient of the 1/e 2 pole, arising at 
two-loop and higher order. 

In our case, the structure of the two Z matrices is 



rkin 




/ z 1 \ 
z 4 z x + z 2 

7 mag 7 mag 7 7 

\Zjrp^ Zj T2 -^l-^mag / 



(14) 



where the renormalization constants Z\, Z% and Z 4 determine the mixing of local, 
dimension-4 heavy-light currents and have been calculated at two-loop order in ||, 
whereas Z mag is the renormalization constant of the chromo-magnetic operator, which 
has been computed at two-loop order in . Here we will calculate the remaining entries 
Zti with the same accuracy. 
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Previous authors have calculated the mixing of the bilocal operators into local oper 
ators at the one-loop order, finding @, |3| 



Zj t1 



O at 



7 mag 
Zj tx 



7 mag 
Z/ T2 



(15) 



Tie Aire 
Using the known expression for the other Z-factors, this information can be used to 
predict the 1/e 2 poles in the two-loop coefficients of the renormalization constants. From 
the second relation in (|13T), we obtain 



7 kin (2) 



ymag (2) 
Zj t1 



7 mag (2) 
J T2 



a s \' 



- 



4:71 J 

4tt/ 
4tt/ 



22 8 
6Cp + -C F C A --C F T F n f 



si 1 



2 ry 

4 



4 
3 



+ — C f Ca 



C F T F n 



4 2 

— C pC 4 — — C F T F n r 

3 3 ^ 



+ 0(« 
+ 0(og), 
+ O(o2). 



(16) 



Here = N , C F = \(N 2 - l)/N and Tp = | are the color factors for an SU{N) gauge 
group, and nf is the number of light-quark flavors. These relations will provide a check 
on our two-loop results. 



3 Two-loop calculation 

To obtain the renormalization constants Zpi at order a 2 sl we calculate the insertions of 
the bilocal operators 0^ n and 0™ ag into the amputated Green function with a heavy 
and a light quark to two-loop order. The relevant diagrams are shown in Figure |l| We 
use dimensional regularization in d = 4 — 2e space-time dimensions and work in the MS 
subtraction scheme. The relation between the renormalized and the bare, regularized 
amputated Green functions is 

T IeQti (v,p, . . .) = (Z h Z q ) 1/2 Z i:j T regtj (v,p, . . .;e) = finite. (17) 

To compute the mixing of the bilocal operators into the local operators with a derivative 
acting on the light-quark field, we only need to keep terms linear in the momentum p of 
the light quark. Because the pole parts are polynomial in the external momenta, we can 
first take a derivative with respect to p and then set p = 0, so that all integrals are of 
propagator type and depend on the single variable uj = v ■ k. However, this method fails 
for some of the diagrams, for which setting p = after differentiation leads to infrared 
(IR) divergences. In these cases, we apply a variant of the so-called R* operation 0, [Tcfl , 
which compensates these IR poles by a recursive construction of counterterms for the 
IR-divergent subgraphs. 

Consider, as an example, the first diagram in Figure |IJ Its contribution is proportional 
to the integral 

D 1 = [d*s dH - ( Jr t \f + P t (S + P \ , • (18) 

J t 2 (t + p) 2 (s + p) 2 (s-t) 2 (v-s + Lu)(v-t + uj) v ; 
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subtract the IR sub divergence of the original diagram, factorizes into an IR counterterm 

dH <ftt+pk _ md/ 2 , ,_ ,w r(d/2)r(d/2-i)r(2-d/2) 
dt t \t + v f ^ pj r(d-i) ' ( j 

and the original diagram with the lines of the IR-sensitive subgraph removed, and with 
t and p set to zero. 

The remaining two-loop tensor integrals are of the general form 

id„Ad+ I W \ / W \ b \i\ ■ ■ ■ b H-n 1 ■■■ L 



d a sd a t 



>vs + uJ \v-t + ujj (-s 2 ) a3 (-t 2 ) Q4 [- (s-t) 2 ] Q5 

= _ vr rf(_2cj) 2( ^ Q3 - a4 - a5)+n+m I£: M ™(v; {a,}) . (21) 
Using the method of integration by parts [JJ |TT[, one obtains the recurrence relation M 
{{d -a x -a*- 2a 5 + n) + a 3 3+ (4" - 5") + a x 1+2"] i^;;^(t;; 

= E^X J l^n( v ;W). (22) 

which can be used to express any two-loop integral in terms of degenerate integrals, 
which have a<i = 0, = or a$ = 0. Here 1 + is an operator raising the index a.\ by 
one unit etc., and \pj] means that this index is omitted. The degenerate integrals can 
be related in a straightforward way to products of one-loop tensor integrals 

Using this technique, we have calculated the pole parts of the two-loop diagrams in 
the 't Hooft-Feynman gauge. The result is a linear combination of terms with the Dirac 
structures corresponding to the local operators 0\ in and 0™2 g defined in (§) and (fT0|). 
The results are summarized in the first two columns of Table |l| In the MS scheme, the 
renormalization scale \i is introduced by the replacement of the bare coupling constant 
with the renormalized one through the relation g^ are = jl e Z g g s with fi = yue 7B,/2 (47r) _1 / 2 . 

The two-loop diagrams in Figure [l] contain sub divergences, which must be subtracted 
by UV counterterms. In addition to the standard one-loop counterterms for the quark 
and gluon propagators and vertices, local operator counterterms are required. To find 
these, we have calculated at one-loop order all insertions of the operators 0^ in and 0™ ag 
into the amputated Green functions with a non-negative degree of divergence. In our 
case, those are the two- and three-point functions with field content qh v , qAh v , h v h v , 
and h v Ah v . We find that in the 't Hooft-Feynman gauge the UV divergences of these 
functions are removed by the counterterms 



C A a s 



Ane 



40 kin - C^ in - 2qT iv D h v - q {%Ip)^Yh v 

+ i J d 4 xT{qTh v (0), h v (iv-D) 2 h v (x)} 

3 3 f 

— -qT v- Ah v — -i J d 4 x T{g Th v (0), h v {iv- D,v- A}h v (x)} 

+ i J d 4 xT{qTh v (0), h v {iD^ A^}h v (x)} 



(23) 
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Table 1: Two-loop and counterterm contributions in units of (a s /47r) 2 



Structure 



Color 



-2uj 



-4e 



-2m 



-2e 



kin 



C f Ca 



4 , ( _ 4 , 8tt 2 \ 1 
c 2 + V 3 + 9 Je 

53 , / 287 _ 4tt 2 \ 1 

6e 2 "T" V 18 3 J e 



_4 

e 

_3 1 

2c 2 c 



--2- 4- 

,2 T 



3c 2 



56 
9c 



97 I 185 
~6c 2 ^ 6c 



16 
3c 2 



32 
3e 



c 2 

C f Ca 



_9_ , /13 , 107r 2 N \ 1 
4c 2 V 3 ~T 9 y e 

_4_ _j_ / 25 5ti£\ 1 

3c 2 ^ V 18 18 J c 



e 



2 

'3c 2 



_1_ 

9c 








7,1 
2e 2 ^ e 



8 
3e 2 

4 
3c 2 



05 



c 2 

C f Ca 
C F T F rif 



_5_ _|_ / 11 , 1 

4c 2 ~T V 6 + 9 je 

_4_ _i_ / 55 \ 1 

3c 2 "r V 18 18 y c 



J_ _ 2 

'c 2 c 



2 

'3c 2 



l_ 

9c 








3 
2e 2 

8 
3c 2 

4 
3c 2 



and 



/'mag 
*"ct. 



47re 



mag + 



mag 
2 



0™ ag + qa^ v YP + a^iv-D h x 



(24) 



Since the two-loop calculation was performed off-shell, some gauge- dependent operators 
must be included in addition to the operators Oj and F , as well as operators that 
vanish by the equations of motion. (The local operators vanishing by the equations of 
motion actually yield no contribution and could be dropped; however, one must keep 
these operators when they appear inside time-ordered products.) The results for the 
counterterm contributions are summarized in the last two columns of Table |l[ When 
these contributions are added to the result for the sum of the two-loop diagrams, all 
nonlocal 1/e divergences proportional to ln(—2uj/{i) and ln(— p 2 //i 2 ) cancel. This is a 
check of our calculation. 

It follows from ( |TTD that the sum of the two-loop diagrams plus their counterterms 
determines the two-loop coefficients of the products — (ZhZg) 1 / 2 Z F %- To obtain the results 
for the renormalization constants Z Fi at order a 2 , we have to account for one-loop 
operator mixing and wave-function renormalization of the external quark fields. Our 
final expressions are 



7 kin 



C F a s 



-4 + 



rre 



4rr 



6 8 



+ T F n f ( - 





8 



^mag 
Zj ti 



C F a s 
Aire 



7 mag 
Z/ T2 



CfCYs 
47TC 



- 1 + 7 1 
An 



C F \ 

Ue 3 



IOtt 2 



Air 



11 2vr 2 



(A 25 Svr 2, 



_ / 4 55 7T 2 
+ Ca 37"18 + 18 



(25) 



They are the main result of this work. Note that the 1/e 2 poles in these expressions 
agree with ([To]). This a nontrivial check of our calculation. 



4 Nontrivial basis transformations 

The results derived in this work are sufficient to calculate, at the two-loop order, the 
operator mixing of bilocal operators into local dimension-4 operators for any choice of 
the Dirac structure V . However, care must be taken when using the results in (GST) for 



a particular choice of T, if instead of the operators in and fllCf ) another set of basis 
operators is employed. Whereas the transformation between one operator basis and 
another is trivial at the one-loop order, it can be subtle at NLO, because in dimensional 
regularization the relations between the operators of different bases may depend on e. 

Consider the general case where the operators {0{\ are expressed in terms of some 
other operators {Qj} by a linear transformation of the form Oi = J2j Rij{ e ) Qj- Depend- 
ing on the choice of T, some of the operators Oi may not be independent, so the new set 
may contain fewer operators than the original one. We define a left-inverse matrix L(e) 
such that L(e)R(e) = 1. It then follows that the matrix Z of renormalization constants 
in the new basis (denoted by a tilde) must be such thatQ 



LZRZ 



= 0. (26) 

poles 



Next, we expand the transformation matrices as R(e) = J2 n Rn(t)£ n and L(e) = 
J2 n L n (e) e n with n > 0, and extract from this result the coefficient of the 1/e pole 
in the matrix Z, which determines the anomalous dimension matrix in the new basis. 
We find 

Z (1) = L Z {1) R Q + [LiZ^Ro + LoZ^Ri 

- 'L^Ro + LqZ^R^LoZ^Ro + . . . , (27) 

where the dots represent terms that do not contribute at two-loop order. 

1 The corresponding relation (19) of || is incorrect; however, this did not affect the results for the 
anomalous dimensions obtained in that paper. 
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Armed with this general result, we now discuss the case of the vector current con- 



sidered in the Introduction. Then the relevant Dirac structures are T = 7" and T 



v 



and the operators Oi are conveniently expressed in terms of the operators Q4...6 defined 

)kir 
T 



in (|). For T = 7 a we have Of n = Q 4 , 0^ n = Q 7 , 0™ ag = Q 9 , and 



01 

0\ 



mag 



(l-e)(l + 2e)Q 4 -4(l-e)Q 5 , 
eQ 4 -2(l-e)Q 5 + (l-2e)Q 6 , 

whereas for T = v a we find Of 1 = Q 5 , 0^ n = Q 8 , 0™ ag = Q w , and 

Or g = -(l-e)(3-2e)Q 5 



(28) 



O 



mag 



-(l-e)Q 5 - 



(29) 



Hence only for the case of insertions of the magnetic operator there appear nontrivial 
basis transformations. 

We thus consider the basis transformation from the six operators consisting of two 
sets of {0™ ag , C>2 iag , 0™ ag } evaluated with T = 7° and T = v a to the new basis spanned 
by the operators {Q4, Q5, Qq, Qg, Qw}- The corresponding transformation matrices -R(e) 
and L(e) can easily be obtained from ( p8|) and (|29|). It is then straightforward to compute 
the matrix Zc corresponding to the operator mixing of Q7...10 into Q4...6, and then using 
the first relation in 
find that 



T3|) to calculate the two-loop anomalous dimension matrix jc- We 
/ 

7c = 

V 



7 kin 





\ 





7i in 





^mag 


7 mag 


mag 
73 





7 mag + 7 mag + ^ag 


J 



(30) 



where 



kin 

7i 



7 mag 



^ag 



mag 
73 



C F a s 

4jT 

4tT 

4ir 

Cp«s 
4tT 



4:71 



2 + ^ 
4tt 



12 + — 

4n 



2 -\ — —— 



C*- - — 



C F — - 



32 


327T 2 


y 


9 


10 


40n 2 


y 


9 



608 16tt 2 



160 



9 



_ /46 IOtt 2 ^ 

My + — 



38 



176tt' 



+ C A 



/236 44tt 2 



4tt 



C F - 



57T 



3 

206 



9 



44^ 
56 

1^ 



9 



+ 



2tT 5 



52 ^ 



This completes the calculation of the anomalous dimension matrix 



(31) 

at two-loop order. 
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5 Wilson coefficients for the vector current 



The Wilson coefficients Bj(fi) in the heavy-quark expansion of the vector current in (|I|) 
obey the renormalization-group equation 



dln/i 



Y ]B{fi) = 0. 



(32) 



where the vector B(fi) contains the ten coefficients Bj(fi). The anomalous dimension 
matrix 7 has been given in (^|), (|7|), (||), and (|30|). Besides the entries of the submatrix j c 
computed in this paper and shown in fl31~l), we need the two-loop anomalous dimension 



7 



mag 



All 



An 



68 „ 52 
— C A -—T F n f 



9 



(33) 



of the chromo-magnetic operator calculated in 
dimensions of local dimension-4 currents, 



as well as the two-loop anomalous 



7i 



72 



73 



CfOs 
All 

CfCVs 
AtT 

CfC^s 

Att 



3 + -^ 

AtT 



3 + 

47T 



2 + 



47T 



-5 + 



13 

3 



"3"/ 

Af_" 
3 / 

9 , 



49 


2tt 2 






If 


+ 3 


41 




3 





10 m 

+ y ^>n/ 



10 



■Tpn/ 



104 2tt 2S 



9 



28^ 



(34) 



which were computed in ||. The coefficient 71 coincides with the universal hybrid 
anomalous dimension of the leading-order heavy-light currents first obtained in 0.0]. 
In addition to the anomalous dimensions, we need the one-loop matching conditions for 
the Wilson coefficients at the scale \i = uiq. In the MS subtraction scheme, they are || 



Ci(m Q ) = Bi(m Q ) 



C F a s 



71 



C 2 {m Q ) = B 2 (m Q ) = -B 3 (m Q ) 



-B 4 {m Q ) = -B 5 {m Q ) = B 6 {m Q ) 



CfOLs 

2tt 

CfOLs 
71 



Cmag(^Q) 



at. 



1 + (C A + C F ) ^ • 

Z7T 



(35) 



The solution of the renormalization-group equation at next-to-leading order is stan- 
dard and described in detail in 0. A subtlety complicating the solution is that the one- 
loop anomalous dimension matrix cannot be diagonalized. We circumvent this problem 
by adding a small diagonal contribution 77 13x3 to the matrix 7^ in (@), taking the limit 
77 — > at the end of the calculation. 
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We quote results for N = 3 colors and n/ = 4 light quark flavors. We first reproduce 
the known NLO results 



d( M ) = B l (/i) = x 



6/25 



1 - 6.243 (0.910 -2A RS ^ aM 



47T 



An 



1 + 8.449 (0.218 + 3A RS ^ ^ 



An 



An 



(36) 



where x = a s (fj,)/a s (mQ). According to @, products of these coefficients determine the 
coefficients -£>7...io- In addition, we obtain the new expressions 



+ x 6/25 lnx 



|_ 9 . 127 ^> + ( 2 . 820 _^A RS 



27 



16 Q.fmo) / 32 A 

_ 3.996 + 0.582 _ Ars 

25 47r V 



+ x 



-3/25 



_±_ 0327 ^q) 

27 4tt 



25 



4 



2.013 — — A R g 



47T 
47T 



47T 



- H + .280^)- 0.993 
9 An 47r 



S 5 (/x) = x 



6/25 



28 , , aAirtn) 

— - + (32.097+ 1.707 In x) s \ Q! 
27 47r 



56 

0.944 - — A RS 



47T 



+ X 



-3/25 



88 aimn) 
— + 11.930 v Q! 
27 An 



30.101 --Ars 



a a (/i) 



4tt 



20 aimn) a s (u) 

— + 0.560 \ g; + 2.458 
9 An An 



BM = x^ 



+ x 



-3/25 



-2 + 7.153 
4 

2.941 

3 An 



An 



1.820 -4A RS ; 



a 8 (fj) 



An 



I 5.246 - -A RS 



An 



+ 



10 



-0.841 e^gl - 1.464 



47T 



47T 



(37) 



All terms proportional to the coupling a s (rriQ) in these expressions, as well as the 
leading-logarithmic rescaling factors, are independent of the renormalization scheme. 
The terms proportional to a s (n) are scheme dependent, however. Here we consider a 
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Figure 2: Different perturbative approximations for the Wilson coefficients in 
the MS scheme: next-to- leading order (solid), leading-logarithmic approximation 
(dashed), and one-loop order (band). 



class of schemes parameterized by the quantity 

A RS = In ^ . (38) 

^RS 

These include all "minimal-subtraction-like" schemes, which are related to each other 
by a change in the renormalization scale. For instance, we have Ams — — 1e + ln47r and 
Ams = by definition, i.e., the above results evaluated with A RS = refer to the MS 
scheme. It follows from ( |32"D that at NLO the scheme-dependent terms in the Wilson 
coefficients are given by 

A%) = ^ 7 r B>). (39) 

In a complete NLO calculation, these terms combine with scheme-dependent terms in 
the hadronic matrix elements of the HQET operators to give a renormalization-group 
invariant answer. 

It is evident from the magnitude of some of the coefficients in (j37|) that the next-to- 
leading corrections are, in some cases, quite large. To illustrate this point, we compare 
our NLO results with the naive one-loop expressions for the Wilson coefficients calculated 
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in j3|, and with the results obtained in leading-logarithmic approximation |2|], where all 
terms proportional to a s /4n are omitted. The different results for the Wilson coefficients 
are shown in Figure |2] as a function of the renormalization scale fi, setting itlq = 4.8 GeV 
corresponding to the 6-quark mass. For the one-loop approximation we show a band 
obtained by varying the coupling between and Q s (mg). We note two important 

observations: (i) The one-loop and leading-logarithmic approximations give rather dif- 
ferent results, which sometimes do not even have the same sign. This shows that the 
one-loop matching corrections at the scale [i = itlq are important, (ii) For low renor- 
malization scales /i ~ 1-2 GeV, some of the coefficients, in particular B^fi) and B 5 (fi), 
are numerically large despite the fact that they vanish at tree level. This shows that 
the effects of running are important, and thus a RG summation of the large logarithms 
a s ln(mg//x) becomes mandatory. The effects of matching and running can only be com- 
bined in a consistent, scheme-independent manner by going beyond the leading order. 
Our NLO results (shown by the solid lines) typically lie between the two approximations, 
except for fi rs itiq, where they agree with the naive one- loop results. 



6 Summary and an application 

We have calculated at two- loop order the mixing of bilocal operators into local dimension- 
4 current operators containing a heavy and a light quark field. When combined with 
results obtained by previous authors, this completes the calculation of the 10 x 10 two- 
loop anomalous dimension matrix governing the renormalization-scale dependence of the 
operators appearing at order I/tuq in the heavy-quark expansion of the weak current 
«7 M (1 — 75)6. We have presented expressions for the Wilson coefficients in this expansion 
at next-to-leading order in renormalization-group improved perturbation theory. We 
find that the next-to-leading corrections are numerically large and must be included for 
a meaningful determination of these coefficients. This makes our results relevant for 
phenomenological applications of the heavy-quark expansion. 

The matrix elements of dimension-4 heavy-light current operators play an important 
role in heavy-flavor phenomenology. They appear, e.g., at order I/tuq in the heavy-quark 
expansion of meson decay constants ||12|| , and of the semileptonic form factors describing 
the exclusive B —>■ ttIv and B —>■ plv transitions ||13|| . The primary motivation for 



our calculation, however, was the recent proposal for extracting the element \V u b\ of the 
quark mixing matrix from the lepton invariant mass spectrum in inclusive B — > X u I v 
decays H . This spectrum can be obtained using a two-step operator product expansion 



(called "hybrid expansion") of a time-ordered product of two heavy-light currents 



Our results will help to reduce the perturbative uncertainty in this expansion, which will 
ultimately lead to an improved accuracy in the determination of \V U b\. 

At NLO in the hybrid expansion, one needs the expression for the combination of 
Wilson coefficients (3I?4 + B e ) derived in the present work. We now present an exact 
result for this combination, evaluating the coefficients with N = 3 colors, but for an 
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arbitrary number of light quark flavors. We find 



1 6 

3s 4 (aO + s 8 (aO = y x2/Po 



where 
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1208 314vr 2 
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28tt 2, 



9 



'941 28tt 2 
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'380 28vr 2 



27 
1 

Po "3 
1 



1 



1177 



9 



27 



Po + 
1 



PI ~ 
-15.476, 

107 



4.098 



214 



2.206 



-6.243. 



(40) 



(41) 



and Po = 11 — |n/ is the first coefficient of the /3- function. The numerical values in 
(p|) refer to n/ = 4. In [IB], we combine this result with an explicit calculation of the 
hadronic matrix elements in the hybrid expansion, finding that the scheme dependence 
parameterized by Ars disappears from the final result for the semileptonic B — > X u I v 
decay rate. 
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